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PERFORMING GEOSCIENCE 
INTERPRETATION WITH SIMULATED 
DATA 

This is a divisional of U.S. application Ser. No. 08/770, 
209, filed Dec. 19, 1996, now U.S. Pat. No. 5,905,657. 

BACKGROUND 

This invention relates to performing geoscience interpre- 
tation with simulated data. 

Geologists, geophysicists, petroleum engineers and pro- 
duction engineers use models, including computerized 
models, of the earth *s shell to plan exploration and produc- 
tion of hydrocarbons and, to a lesser extent, other minerals. 
As hydrocarbons become more and more scarce, the accu- 
racy of the computerized models becomes increasingly 
important to limiting the cost of locating and producing 
hydrocarbons and the associated cost of hydrocarbon 
products, such as gasoline and heating oil. 

Existing modeling systems use an "inversion" technique 
to model geoscience structures from acquired data. Acquired 
data, such as velocity measurements or data collected 
through seismic instrumentation, is "inverted" and used to 
predict the location of subsurface structures and the physical 
properties of those subsurface structures. 

Once the inversion process produces predictions of the 
locations and physical properties of the subsurface 
structures, they are modeled geometrically. Such models 
must be capable of representing complicated geometrical 
shapes, such as the shape of a salt dome or salt wall or the 
shape of a group of folded sedimentary beds. 

Commonly, modeling systems use a grid system to rep- 
resent subsurface structures. In a grid system, the subsurface 
structures are represented by points in space regularly 
spaced in the x-y plane and having a single value in the z 
plane. The values corresponding to those subsurface struc- 
tures between the regularly-spaced points are determined by 
interpolation. In a model for a large subsurface region, the 
number of points in space can be numerous, requiring a great 
deal of storage, even if the subsurface structures are uniform 
in topography and physical properties across that space. 

Another modeling system that has been applied to the 
problem of modeling subsurface structures is constructional 
solid geometry ("CSG"). Under this approach, two or more 
solid primitives (such as cylinders or spheres) are combined 
to form a more complicated solid. The process begins with 
objects that are known to be solid and adds or subtracts other 
solids to create new shapes. For example, subtracting a 
cylinder from a box will produce a box with a hole in it. 

SUMMARY 

In general, in one aspect, the invention features a method 
for analyzing geological data stored in a geoscience model 
on a magnetic media, comprising building a simulation input 
model from the geoscience model. 

Implementations of the invention may include one or 
more of the following. The method may further comprise 
acquiring data; interpreting the acquired data to produce the 
geoscience model; applying a simulator to the simulation 
input model to produce synthetic data; comparing the 
acquired data to the synthetic data to produce a difference; 
and editing the geoscience model to reduce the difference. 
The simulator may comprise a simulation program used to 
develop acquisition equipment. Building the simulation 
model may comprise producing a finite element mesh. 
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Building the simulation model may comprise producing a 
finite difference grid. Building the simulation input model 
may comprise producing a tessellated model. The simulation 
input model may have a different format than the geoscience 
5 model. 

In general, in another aspect, the invention features a 
method for analyzing geological data representing a subsur- 
face region, the geological data being stored in a geoscience 
model on a magnetic media, comprising constructing a 
10 boundary for dividing the region into a first sub-region and 
a second sub-region; storing on the magnetic media the 
shape of the boundary as a parametric function whose 
parameter density can vary. 

In general, in another aspect, the invention features a 
15 method for analyzing geological data sampled from a sub- 
surface region, the geological data and an analysis of the 
geological data being stored in a geoscience model on a 
magnetic media, the method comprising dividing the region 
into sub-regions, in each of which a material property varies 
20 without discontinuities. 

Implementation of the invention may include one or more 
of the following. The method may further comprise describ- 
ing the variation of the material property with a parametric 
function. 

25 

In general, in another aspect, the invention features a 
computer system for analyzing geological data stored in a 
geoscience model on a computer-readable magnetic media, 
comprising means for building a simulation input model 

3Q from the geoscience model. 

Implementations of the invention may include one or 
more of the following. The computer system may further 
comprise means for acquiring data; means for interpreting 
the acquired data to produce the geoscience model; means 

35 for applying a simulator to the simulation input model to 
produce synthetic data; means for comparing the acquired 
data to the synthetic data to produce a difference; and means 
for editing the geoscience model to reduce the difference. 
The simulator may comprise a simulation program used to 

40 develop acquisition equipment. The means for building the 
simulation input model may comprise a means for producing 
a finite element mesh. The means for building the simulation 
data input model may comprise a means for producing a 
finite difference grid. The means for building the simulation 

45 input model may comprise a means for producing a tessel- 
lated model. The simulation input model may have a dif- 
ferent format than the geoscience model. 

In general, in another aspect, the invention features a 
computer system for analyzing geological data sampled 

50 from a subsurface region, the geological data and an analysis 
of the geological data being stored in a geoscience model on 
a computer-readable magnetic media, comprising means for 
constructing a boundary for dividing the region into a first 
sub-region and a second sub-region; means for storing on 

55 the magnetic media the shape of the boundary as a para- 
metric function whose parameter density can vary. 

In general, in another aspect, the invention features a 
computer system for analyzing geological data sampled 
fr6m a subsurface region, the geological data and an analysis 

60 of the geological data being stored in a geoscience model on 
a magnetic media, the computer system comprising means 
for dividing the region into sub-regions, in each of which a 
material property varies without discontinuities. 

In general, in another aspect, the invention features a 

65 computer program, residing on a computer-readable 
medium, for analyzing geological data sampled from a 
subsurface region and an analysis of the geological data 
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stored in a geoscience model on a magnetic media, com- The geoscience model 12 is a permanent repository of the 

prising instructions for causing a computer to build a simu- interpretation which is stored in a general format so that 

lation input model from the geoscience model. input for a variety of simulators can be generated from it. 

Implementations of the invention may include one or The geoscience model is editable so that interpreters can 
more of the following. The computer program may further 5 quickly converge the model to a consistent interpretation, 
comprise instructions for causing the computer to acquire Once the geoscience model exists and has been verified as 
data; interpret the acquired data to produce the geoscience consistent with existing data, it can be used for a variety of 
model; apply a simulator to the simulation input model to predictive and planning purposes such as survey planning, 
produce synthetic data; compare the acquired data to the estimating reservoir performance, and well completion plan- 
synthetic data to produce a difference; and edit the geo- , 0 ning. 

science model to reduce the difference. The simulator may The same technologies needed to support simulation 

comprise a simulation program used to develop acquisition based interpretation can be used to improve the process for 

equipment. The instructions which cause the computer to designing acquisition equipment. Equipment models placed 

build the simulation input model may comprise instructions within geoscience models may be used to build simulation 

which cause the computer to produce a finite element mesh. input models for simulators. Simulation outputs are ana- 

The instructions which cause the computer to build the 15 lyzed and compared to design specifications. Design alter- 

simulation input model may comprise instructions which natives and parameter variations can be quickly tested and 

cause the computer to produce a finite difference grid. The optimized prior to building the first prototype, 

instructions which cause the computer to build the Simula- ^ same technologies needed to support simulation 

tion input model may comprise instructions which cause the based inter p r e ta tion can be used to improve the quality of 

computer to produce a tessellated model. The simulation 20 ed process of survey planning. 

input model may have a different format than the geoscience ^ a ■ ■ ? • ai 1 /■ 

m ^ de j J to Source and receiving equipment models placed in geo- 

T " . . . . . - science models may be used to generate simulation input 

In general, in another aspect, the invention features a _ j 1 r - w c . , ^ _ ^ M , , , 

computer program, residing on a computer-readable m^els for simulators. Smulatton outpute may be analyzed 

medium, for analyzing geological data saVnpled from a 25 and used to °P tunize * e P Iacement of sources ™ d 

subsurface region, the geological data and an analysis of the receivers. 

geological data being stored in a geoscience model on a ° ne P** of geoscience model is a geometric repre- 

magnetic media, comprising instructions for causing a com- sentation of a portion of the earth's subsurface. In geo- 

puter to construct a boundary for dividing the region into a science applications, known geometric primitives, such as 

first sub-region and a second sub-region; store on the 30 s P heres and c y linders * « less usefuI than * other a PP uca - 

magnetic media the shape of the boundary as a parametric tions - underground structures may have started out 

function whose parameter density can vary. ±eir geological life capable of being described as simple 

In general, in another aspect, the invention features a geometric primitives, over time :most have become distorted 

computer program, stored on a computer-readable medium, ™ d m not ^ ^presented using a grid system or a 

for analyzing geological data sampled from a subsurface 35 constructional solid geometry system. For example, if a set 

region, the geological data and an analysis of the geological of sedimentary beds folds over on top of itself a single x-y 

data being stored in a geoscience model on a magnetic Pomtwill have two or more z values for a given bed. Further, 

media, the computer program comprising instructions for the sh *P e * « vei 7 com P hcated ^ not e ^ represented 

causing a computer to divide the region into sub-regions, in usin g primitives. 

each of which a material property varies without disconti- 40 ^ role of geoscience model is to represent the shape 

nu j t j es and properties of a system under study so that simulation 

programs may be used to enhance understanding of the 

DESCRIPTION system, as shown in FIG. 2. A variety of 3D modeling 

FIGS. 1-3 are block diagrams. applications 20 including velocity modeling 22, borehole 

FIGS. 4a-c y 5a-c and 6a-g are representations of mod- 45 interpretation 24, seismic planning 26, and others 28, share 

eled items. a common geoscience model 12. A simulation input model 

Simulation based interpretation ("SBI"), illustrated in generator 30 generates from the geoscience model inputs 

FIG. 1, is an alternative to the inversion process for pro- required for a variety of simulators 32, such as migration 34, 

ducing consistent models from acquired data. SBI begins by sonic 36 > nuclear 38, electromagnetic 40 and other types of 

gathering acquired data 10 from the geologic structures 50 simulators 42. 

being modeled. The acquired data is interpreted and the A simulation input model uses the data format required by 

resulting interpretation is captured as a geoscience model a simulation program. Each simulation program has its own 

12. The model is then tested by simulating an original specialized input data format. These formats are specified 

acquisition experiment to produce synthetic data 14. For for a variety of reasons including, but not limited to, run time 

example, the simulation programs used to develop acquisi- 55 efficiency, data storage efficiency, numerical methodology, 

tion equipment can be used to produce the synthetic data. and design history. 

Differences between the synthetic data 14 and the acquired The current generation of simulation programs are based 
data 10 are determined and used to adjust the model. primarily upon the finite element, the finite difference, and 
Repeated iterations of these steps may cause the model to ray tracing methods. The simulation input model generator 
converge toward a consistent model in which the differences 60 provides automated procedures for generating finite differ- 
between the synthetic and the acquired data are reduced. ence grids, finite element meshes and tessellated models for 
Making a single interpretation consistent with several mea- ray tracing that may be specialized to a particular simulation 
surements from a variety of experiments at one time con- program's requirements. 

tinues to improve the quality of the interpretation. Trained The domain of the geoscience model is defined by the 

interpreters can use such a system interactively to build 65 simulation programs it is required to support. Despite the 

interpretations which can be relied upon for prediction many different physical phenomena that the simulators 

purposes. model, they all have the following in common: 
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(a) the continuum assumption; 

(b) a representation of distributed material properties; 

(c) a set of shared simulation techniques including the 
finite element, finite difference, and ray tracing meth- 
ods. 

The continuum assumption states that at a sufficiently 
large scale material properties, such as porosity and velocity, 
of a sub-region, which is a portion of a region, can be 
approximated by macroscopic properties that vary smoothly 
throughout the sub-region. The locations of discontinuities 
in the material properties, such as at a fault, are called 
"boundaries." Material property fields and their discontinui- 
ties can be represented in a geoscience model which con- 
tains the following information: 

(a) Shape, which is defined to be the location of the 
discontinuity boundaries; 

(b) Topology, which is defined to be the connections 
between boundaries and sub-regions; 

(c) Material properties, which are defined to be smoothly 
varying values within and on the boundaries of each 
sub-region and boundary object; and 

(d) Simulation input models, which are defined to be an 
extraction of the information in a geoscience model in 
a form suitable for use by a simulation package. 

The geoscience model and the simulation input models 
are distinguished so that any number of simulation packages 
can be supported by a single geoscience model. The simu- 
lation input models support simulation based interpretation 
applications. 

Irregular space partition ("ISP") models are a represen- 
tation of smoothly varying material property fields which 
may contain internal discontinuities that occur along distinct 
boundaries. A discontinuity in a material property occurs 
where the material property value differs by a finite amount 
on either side of a discontinuity boundary. The location of 
the discontinuity boundaries make up the shape of the 
model. ISP models explicitly represent the shape of the 
discontinuity boundaries. These boundaries may lie any- 
where in an infinite space, have any shape, and intersect one 
another to partition space into sets of distinct sub-regions. 
Within any one sub-region, the material property can be 
represented by a smoothly varying function. 

The boundaries can be highly irregular in shape and can 
comprise portions where their shape is fairly uniform (such 
as a plane or a surface that approximates a plane) and other 
portions where the boundary is changing rapidly (such as 
where folding occurs). The uniform portions of the bound- 
aries can be represented accurately with fewer parameters. 
For example, a plane can be represented by three points. In 
the portions where the boundary is changing rapidly, 
however, it is necessary to use more parameters to achieve 
an accurate representation. 

For example, a "hybrid" surface may be made from both 
a grid and a mesh. The grid represents the bulk of the surface 
(typically the interior) where a high density of points is not 
necessary. The mesh represents the surface boundary areas 
and surface intersections with other surfaces, curves or 
points where high density is required. Thus, the surface is 
represented with a high density of parameters only where 
such a representation is necessary. This allows the system to 
store large surfaces efficiently while maintaining higher 
accuracy of surface boundaries, when compared to regular 
grid representations. 

ISP models are also used to model manufactured objects 
such as acquisition, drilling, and truck equipment, which are 
represented as solid models. Points in a solid model are 



classified as being either in or out of the modeled object. In 
such cases, boundaries are not discontinuities between two 
regions but the demarcation of the end of the model. In a 
sense, there is nothing on the other side of a boundary. Solid 

5 models do not contain internal boundaries. 

A three-dimensional simulation-based application 50 
includes three-dimensional modeling 52 which comprises 
four basic geoscience model building steps divided into 
geometry modeling 54 and simulation input model genera- 

10 tion 56, as shown in FIG. 3. The sequence of shape editing 
58, topology model building 60, material property editing 
62, and simulation input model generation 56 (which com- 
prises finite element meshing/finite difference gridding/ 
tessellation for ray tracing) defines the scope of geoscience 

15 modeling. 

As discussed above, a geoscience model may consist of a 
combination of a large number of simple geometry primi- 
tives. Each geometry primitive is initialized, shaped, and 
positioned interactively by a user in relationship to the rest 

20 of the geoscience model. These geometry primitives include 
such shapes as lines, circles, planes, free-form surfaces, and 
solid sub-volumes. 

The shape of a geometry primitive is defined to be the set 
of points it occupies in space. For example, a circle of radius 

25 R, lying in a plane parallel to the x-y plane and centered on 
the point [a, b, c] is the set of points given by the function 



30 



x=a+R cos (theta) 
y=tH-R sin (theta) 



where theta is in the range 0<=theta<2 n. 
These equations define a map from theta space to [x, y, z] 

35 space; that is, for every value of theta there is a unique value 
of [x, y, z]. Since theta can be any of an infinite number of 
points the number of points in the circle is also infinite. 

The values {a, b, c, R} are called the parameters of the 
circle. Together, they completely specify all the points that 

40 define the circle. The functions which represent geometry 
primitive shapes are all parametric, single valued, continu- 
ous functions. These functions have well behaved math- 
ematical properties upon which the algorithms used to build, 
edit, and query complicated geoscience models depend. 

45 Parameterization allows an infinite number of points to be 
specified by a discrete number of values. Interfaces can be 
built that conveniendy allow the users to edit the parameter 
values. For example, a pick and drag window-based inter- 
face allows a user to select the circle's center point and drag 

50 it through space. The code that implements this behavior has 
only to convert the mouse's position into a spatial (x, y, z) 
location and use those values to change the shape of the 
circle, e.g. a=x, b=y, c=z. 

Shape editing is the interactive process of modifying a 

55 shape's parameters until it is in a form suitable for use with 
the rest of the geoscience model. In the circle example, the 
object has only 4 parameters. In general, an object may have 
many more parameters. For example, a free-form surface 
represented as a triangulated mesh can have thousands of 

60 parameters. The application which helps a user edit the 
shape of an object with a large number of parameters may 
very well be a complicated application in itself. But all shape 
editing applications are dedicated to creating new geometry 
primitives and setting their parameter values. 

65 There are two kinds of shape editing applications: those 
that extract shapes from data, and those that build shapes 
from scratch. Extracting shape from data is a user directed 
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semi-automated process, while building shape definitions 
from scratch is interactive. 

ISP and solid models are built by combining sets of 
geometry primitives and recording where the primitives 
intersect one another. In three-dimensional ISP building, a 
user creates a number of surface boundary shapes and uses 
them to sub-divide space locally to define a set of sub- 
volumes in a single model. Where boundaries intersect one 
another, new intersection geometry objects are created to 
represent the intersection shape and to record the connection 
among the boundaries. When sets of boundaries interact to 
separate one region of space from another, a sub-region 
object is created for each distinct sub-region. Groups of 
sub-regions may be used to represent application features 
such as formations and invasion zones, and boundaries may 
be used to represent application features like horizons, 
faults, and borehole boundaries. Not every object in a 
geoscience model is used to represent an application feature, 
but every application feature with a geometric description is 
associated with a geometry object. Geometry objects repre- 
senting the shape of an application feature are called geom- 
etry features. 

Building solid models is very similar to building ISP 
models. In solid modeling, the user builds a sequence of 
solid shapes and combines them together with the standard 
boolean operators of union, subtract, and intersect to create 
accurate representations of complicated machined metal 
solids. The side effect of editing the topology structure 
representation when combining primitives into ISP and solid 
models is automated. Sequences for building ISP and solid 
models may use the subdivide, insert, and subtract operators, 
which are just a few of the topology editing operators to be 
found in a general geometry modeling application and which 
are described in detail in co-pending application Ser. No. 
08/772,082, U.S. Pat. No. 6,128,577 entitled "Modeling 
Geological Structures and Properties," incorporated by ref- 
erence. 

An example of ISP modeling is shown in FIGS. 4a-c. A 
volume 70 may be subdivided by two surfaces 72 and 74, as 
shown in FIG. 4a. Another volume 76, representing, for 
example, a salt dome, may be inserted into the subdivided 
volume 78 as shown in FIG. 4b. The result is a subdivided 
volume 80 containing a variety of complicated shapes, as 
shown in FIG. 4c. An analyst can mark interesting geometry 
objects such as the horizons, formed by surfaces 72 and 74, 
or the salt body 76. 

An example of solid model building is illustrated FIGS. 
5a through 5c. A series of sub- volumes 82a-d is subtracted 
from an initial shape 84, as shown in FIG. 5a to produce the 
volume 86, illustrated in FIG. 5b. Additional sub-volumes 
SSa-c are subtracted from the volume, as shown in FIG. 5B, 
resulting in a sub-volume 90, shown in FIG. 5C. 

Material property editing is the process of annotating a 
geoscience model with material properties. Material prop- 
erty fields, like shape, are represented by functions. A 
parametric function is used to map each point in a modeling 
image space to a material property value. Any continuous 
function may be used. Some of the more common functions 
include: constants, constant gradients, and free-form func- 
tions (such as triangulated surfaces, poly-lines, NURBS, 
bi-cubic patches, etc.). 

The shape and topology information in a geoscience 
model is used to model material property discontinuities. 
Each material property discontinuity is marked by a geom- 
etry boundary creating sets of sub-regions in which the 
material property field is smoothly varying. Each sub-region 
is associated with its own function for representing a com- 
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mon material property field. Thus, material property discon- 
tinuities are represented by using different material functions 
on either side of a discontinuity boundary. It takes two steps 
to evaluate the material property at a point within a geo- 
science model: 

(a) find the sub-region containing the point and its asso- 
ciated material property function; and 

(b) evaluate the sub-region's material property function at 
the required point. Any number of material property 
fields can be represented in a geoscience model by 
storing one parametric function for each material prop- 
erty in each sub-region. 

Material property fields may also be defined on region 
boundaries. These are known as "boundary material prop- 
erties" as opposed to "region material properties." Boundary 
material properties can be associated with the entire bound- 
ary or with just one side of the boundary. Oriented boundary 
material properties may be used to describe boundary value 
problems. 

The value of a material property at a point is a tensor. The 
tensor may be represented by either a scalar, vector, or 
matrix quantity. It is a description of properties with intrinsic 
geometric properties, e.g., it remains invariant under coor- 
dinate transformations. The components of a tensor are 
represented in a chosen coordinate system. Changing the 
coordinate system changes the values of the coordinates. 
Relationships between tensor quantities as defined by fun- 
damental physical laws are independent of coordinate trans- 
formations. For example, Newton's law defines the relation- 
ship between acceleration, A, and force, F, as, F=mA, which 
defines the quantity mass, m. This relationship is true 
independent of which coordinate system is chosen to 
describe the components of the tensor quantities F, m, and A. 
In this system, there is no special accommodation for 
tensors. The components of a high order tensor are stored as 
a set of scalar values and the applications manage the 
constraints between component values. 

Material property editing, like shape editing, occurs in 
two ways: building material property fields from data and 
from scratch. Material property editing consists of: 

(a) Specifying the set of material properties to be repre- 
sented in the geoscience model. 

(b) Selecting and assigning a parametric function for each 
material property 

in each sub-region, and optionally 

in each boundary, and 

on each side of a boundary. 

(c) Editing the parameters of the material property para- 
metric functions as desired. 

Representing a material property field in an ISP geo- 
science model is done by using functions. Functions are used 
to represent the boundary shapes which mark the disconti- 
nuities within the material property field and map points 
from a parameter space to a shared image space. Each 
boundary has its own parameter space but all boundaries 
share a common image space. The parameter space is the 
domain of the shape function and the shared image space is 
its range. The shared image space is used as the domain for 
a function which represents a material property field. Every 
point in the shared image space is projected into a material 
property value by the material property function. 

The steps used to subdivide a volume and assign material 
properties to the volume are illustrated in FIGS. 6a-g. A first 
horizon 100 is represented in parameter space, as shown in 
FIG. 6a. A function, f, 102, provides a mapping from 
parameter space to region space. The representation 104 of 



